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THREE-DIMENSIONAL SIMULATIONS OF MAGNETOHYDRODYNAMIC WAVES IN MAGNETIZED SOLAR 

ATMOSPHERE 

Vigeesh, G. 1,2 ; Fedun, V. 3 ; Hasan, S. S. 2 ; Erdelyi, R. 3 
Abstract 

We present results of three-dimensional numerical simulations of magnetohydrodynamic wave prop- 
agation in a solar magnetic flux tube. Our study aims at understanding the properties of a range of 
MHD wave modes generated by different photospheric motions. We consider two scenarios observed 
in the lower solar photosphere, namely, granular buffeting and vortex-like motion, among the simplest 
mechanism for the generation of waves within a strong, localized magnetic flux concentration. We 
show that, granular buffeting is likely to generate stronger slow and fast magneto- acoustic waves as 
compared to swirly motions. Correspondingly, the energy flux transported differ as a result of the 
driving motions. We also demonstrate that the waves generated by granular buffeting are likely to 
manifest in stronger emission in the chromospheric network. We argue that different mechanisms 
of wave generation are active during the evolution of a magnetic element in the intergranular lane, 
resulting in temporally varying emission at chromospheric heights. 

Subject headings: Magnetic fields, Magnetohydrodynamics (MHD), Waves, Sun: atmosphere, Sun: 
surface magnetism 



1. INTRODUCTION 

The highly dynamic solar photosphere is a rich source 
of waves, which are believed to be one of the plau- 
sible candidates responsible for heating the upper so- 
lar atmosphere (see a recent review by Taroyan & 
Erdelyi 2009, in this context). Magnetic fields and 
magnetohydrodynamic (MHD) waves play an important 
role in the solar atmosphere's energy budget. With 
recent high-resolution observations using both space- 
borne (Hinode/SOT, SoHO, SUNRISE) and ground- 
based (SST, DST/ROSA) telescopes and simulations by 
various groups (e.g. Wedemeyer-Bohm & Rouppe van der 
Voort 2009; Shelyag et al. 2011b) showing that the solar 
photosphere harbours vortex motions at the downdraft 
regions of convective cells, here we aim to investigate the 
efficiency of these motions in the generation of a range of 
MHD modes and their transport of energy. Using three- 
dimensional simulations, we show that the excited MHD 
modes and the wave energy transported by these modes 
vary with different excitation scenarios. 

Several authors have carried out studies of wave prop- 
agation in magnetic flux tubes. Some of the recent work 
on magnetoacoustic wave propagation by e.g. Hasan & 
van Ballegooijen (2008), Khomenko et al. (2008), Vigeesh 
et al. (2009), and Vigeesh et al. (2011) considered mag- 
netic network fields with a range of typical properties. 
They have highlighted the importance of mode coupling 
in the transport of energy by magneto-acoustic waves in 
a magnetised atmosphere. These studies were carried out 
in 2D, and they do provide a reliable picture of the actual 
processes taking place in the real atmosphere. A more 
recent work by Fedun et al. (2011a) considers an open 
flux tube in the two-dimensional solar atmosphere with 
combination of realistic temperature profiles correspond- 
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ing to different height ranges. These investigations bring 
out the importance of the transition region in influencing 
the propagation of magneto-acoustic waves. While a lot 
of information can still be unveiled in 2D with realistic 
atmosphere, similar to work carried out by Fedun et al. 
(2011a), a 3D approach to the problem is more appealing 
in terms of physical processes, although more expensive 
in terms of computational resources. One of the major 
new features by allowing an extra dimension to the prob- 
lem is the possibility of the excitation of the "intermedi- 
ate" Alfven wave. The propagation of MHD waves in a 
three-dimensional atmosphere by the driving motions at 
the lower boundary has been studied by Malins & Erdelyi 
(2007), Erdelyi et al. (2007), Fedun et al. (2009). They 
considered a uniform background magnetic field embed- 
ded in a realistic hydrodynamic model of a quiet Sun and 
have shown that magnetic field channels the photospheric 
disturbances all the way into the solar corona. These 
studies did not consider Alfven waves due to the non- 
torsional nature of the driving motion. Further studies 
of MHD simulations in three-dimensional models (Fedun 
et al. 2011b) have, however, shown that the vortex mo- 
tions in the solar photosphere (Shelyag et al. 2011a) do 
generate different types of MHD modes, including tor- 
sional Alfven waves. In a recent paper, Fedun et al. 
(2011c), have demonstrated that torsional Alfven waves 
can reveal the magnetic structuring of chromospheric fea- 
tures, since the magnetic flux tubes act as a spatial fre- 
quency filter (Jess et al. 2009; Verth et al. 2010). These 
studies (e.g. Vigeesh et al. 2011; Fedun et al. 2011c) high- 
light that the various MHD modes are a valuable tool 
for the magneto-seismic studies of the lower solar at- 
mosphere. A considerable leap forward connecting the 
swirly lower solar atmosphere to the corona is made by 
Wedemeyer-Bohm et al. (2012) where the role of Alfven 
waves in channelling photospheric energy to the corona 
is demonstrated both observationally and confirmed by 
modelling. 

High-resolution observations using the Swedish Solar 
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Telescope (SST) (Bonet et al. 2008) have revealed con- 
vectively driven vortex flows in the photosphere of the 
Sun. These flows result from downdrafts that occur in 
the regions where the flows converge at the boundary 
of convective cells. Bonet et al. (2008) found that the 
associated magnetic concentrations, visible as magnetic 
bright points (MBP) follow a spiral pattern as they are 
carried along by the vortical flow. Wedemeyer-Bohm 
& Rouppe van der Voort (2009) have detected simi- 
lar features in the core of the chromospheric line Ca II 
(854.2 nm) as a likely result of the disturbance caused 
in the magnetic fields due to the twisting of the under- 
lying photospheric magnetic field. New results from the 
SUNRISE experiment indicate that the vortex flows are a 
very common phenomena in the solar atmosphere (Bonet 
et al. 2010). They are formed in localized convective 
downdrafts as a result of angular momentum conserva- 
tion. These observations also reveal vortex tubes (Steiner 
et al. 2010) with vorticity directed horizontal and are 
considered to cause asymmetric brightenings in the gran- 
ular boundaries. The direct detection of torsional Alfven 
oscillations associated with MBPs, as reported by Jess 
et al. (2009), can be considered a likely consequence of 
the swirling of the magnetic-flux tubes as a result of the 
vortical motions. 

Vortex-like features are also seen in large-eddy simula- 
tions of Shelyag et al. (2011b) pointing to the fact that 
these vortical motions are common feature in the magne- 
tised photosphere as a result of a convective driver. This 
numerical prediction is now observationally confirmed by 
Wedemeyer-Bohm et al. (2012). Using realistic three- 
dimensional radiative hydrodynamic simulations, Kiti- 
ashvili et al. (2011) show that interactions of vortices 
in the intergranular lanes excite acoustic waves into the 
overlying layers. Moll et al. (2011) have also reported the 
occurrence of vortices in surface convection simulations. 

In this paper, we model individual flux tubes located 
in inter-granular lanes in an atmosphere representing the 
empirical temperature profile of the solar chromosphere. 
Convective motions are thought to excite waves in this 
medium, modelled here by means of a driver located in 
the photosphere. Two different driving mechanisms are 
considered, viz. transverse uni-directional and torsional 
motions. The driving generates various MHD modes 
within the flux tube and acoustic waves in the ambient 
medium. We study the propagation and other properties 
of MHD modes generated by these drivers. 

Section 2 discusses the construction of the initial 
magneto-hydrostatic equilibrium model and the proper- 
ties of the model. The excitation mechanisms are ex- 
plained in Sec. 3. The numerical methods and boundary 
conditions are described in Sec. 4. In Sec. 5 we discuss 
the results and the conclusions are presented in Sec. 6. 

2. INITIAL MAGNETO-HYDROSTATIC MODEL 

We model a strong, axially symmetric magnetic flux 
tube in a stratified solar model atmosphere whose lower 
boundary is located at photospheric levels. The horizon- 
tal distribution of the magnetic field components (B x and 
By) across the flux tube at different leveles are shown in 
Fig. 1. The initial magnetic field strength on the axis at 
z = km is 1435 G, which decreases to 161 G at the top 
of the box almost to a value close to the magnetic field 
strength in the ambient medium. The Alfven speed (va) 
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Fig. 1. — B x and B z components of the magnetic field as a func- 
tion of the horizontal distance at the following heights: z = km 
(solid curve), z = 0.5 Mm (dotted curve), and z = 1 Mm (dashed 
curve) of flux tube. 

in the ambient medium changes by four orders of mag- 
nitude as we ascend from z = Mm to z = 1 Mm, but 
it changes only from 10.9 km s _1 to 77.7 km s _1 inside 
the magnetic flux tube. The initial model is constructed 
in a physical domain of size: 1 Mm x 1 Mm x 1 Mm and 
has 100 x 100 x 100 grid points, respectively. The model 
considered in our paper is an isolated flux tube which 
in a Cartesian frame is symmetric with respect to the 
x = Mm and y = Mm planes. The construction of 
the model is based on a 3D generalisation of the method 
that was presented in Hasan et al. (2005). Close to the 
base {z = Mm) the flux tube has a nearly circular cross- 
section with radius smaller than the half-width of the 
computational box, L = 0.5 Mm. For low heights, the 
flux tube is surrounded by a nearly field-free medium. It 
expands with height and its radius becomes comparable 
to L. To simulate the effect of neighbouring flux tubes, 
we assume that the normal component of the magnetic 
field vanishes at the side walls of the computational do- 
main. More details on the 3-D MHD equilibrium model 
are given in the Appendix. The characteristic physical 
parameters of the magnetic flux tube are shown in the 
Table 1. The values in brackets correspond to the am- 
bient medium. The temperature of the magnetic flux 
tube increases with height, correspondingly increasing 
the sound speed (c$) from 7.7 km s _1 to 8.8 km s -1 . 
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TABLE 1 

Equilibrium model characteristics for the flux tube. 
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1.6 
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77.7 
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x 10" 7 ) 
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(39.2) 


(141) 


(0.06) 


z=0 Mm 


4768 


1.3 


X 10" 4 


4.2 x 10 3 


7.1 


10.9 


1435 


0.5 


(4766) 


(4.0 


x 10- 4 ) 


(1.2 x 10 4 ) 


(7.1) 


(0.003) 


(0.77) 


(5.1 x 10 6 ) 



To aide the reader in visualizing the model, we plot 
a three-dimensional representation of the magnetic flux 
tube in Fig. 2 by visualizing field lines that pass through 
points at distance of r = 0.35 Mm from the axis and 
height of z = 1 Mm. The light gray lines that mark 
the field lines can be thought as curves that form the 
magnetic flux iso-surface, since the flux tube that we 
construct here is axisymmetric. In the figure we also 
highlight a specific field line (thick dashed curve) and 
the unit vectors corresponding to different directions at 
an arbitrary point P on this field line. Further details will 
be discussed in Sect. 5 where we introduce the various 
components of velocity with respect to the flux tube. In 
order to have a clear picture of the shape of the tube 
and the surface of equipartition between thermal energy 
density and magnetic energy density, we plot in Figure 3 
the magnetic field lines of the flux tube on an x — z plane 
at y = Mm. The thick curve traces the plasma (3 = 1 
surface (/? is the ratio of gas to magnetic pressure) which 
represents the equipartition level. When viewed in three 
dimensions, the angle between the magnetic vector and 
the (3 = 1 surface normal varies between 0° and 90° in 
the entire computational domain. 

In Fig. 4 we plot the inclination angle of the magnetic 
field vector and the (3 = 1 surface normal at three differ- 
ent x — z planes corresponding to y = Mm, 0.2 Mm, and 
0.4 Mm. The temperature profile within the flux tube 
and in the ambient medium of the equilibrium model is 
shown in Fig. 5. 

A typical magnetic flux tube model like the one de- 
scribed above can support the well known three MHD 
modes. They are the slow, fast magneto-acoustic and the 
intermediate modes, named so according to the speed of 
propagation respective to the local Alfven speed (va)- 
In a low-/? medium, the slow wave has a phase speed 
less than or equal to the sound speed, since the wave 
is essentially driven as a result of gas pressure perturba- 
tions. Hodographs or Friedrich diagrams for MHD waves 
show that the phase speed of the slow wave equals to the 
sound speed (low-/?) for a direction along the field line 
and drops to zero as a function of the angle between the 
propagation vector and the magnetic field vector. On 
the contrary, the fast wave propagates independent to 
the field line with a maximum speed equal to vp, the 
fast speed. 

Figure 6 shows the Alfven speed (va), sound speed (c$) 
and the fast speed (vp) defined by: 

v F = \]c 2 s + v% (1) 
and the tube speed (cp) given by: 

/ °S 2v A 2 , v 

V c s 2 + v A 2 
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Fig. 2. — Three-dimensional representation of the magnetic 
flux tube. The magnetic flux iso-surface is marked by the field 
lines (thin gray lines) passing through points at a distance of 
r =0.35 Mm from the axis at a height of z =1 Mm. The thick 
dashed line represents a single field line on this iso-surface. 
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Fig. 3. — Two-dimensional representation of the flux-tube. The 
thin lines mark the magnetic field lines on a x—z plane at y =0 Mm. 
The field lines that trace the iso-surface shown in Fig. 2 is shown 
as thick lines. The (3 =1 contour is shown as thick dashed line. 

as a function of height on the axis of the flux tube, on 
the representative field line and in the ambient medium 
of the model. This plot helps us to identify the differ- 
ent regimes of magneto-acoustic wave propagation in the 
model. To clearly understand the wave behavior at dif- 
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Fig. 4. — The inclination of the (3 =1 surface normal to the local 
magnetic field vector. The thick curve is for a plane at y =0 Mm. 
The dotted line is for a plane at y =0.2 Mm and the dashed line is 
for y =0.4 Mm. 
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Fig. 5. — Temparature as a function of height. The thick curve 
represents is for the axis and the dotted curve is for the ambient 
medium. 



ferent locations in the domain, we separate it into six 
different regions. In the following, represents the an- 
gle between the propagation direction and magnetic field 
vector. 

On axis: The light-grey shaded region (Region I) corre- 
sponds to the propagation of fast wave on the axis of 
the flux tube. The upper limit of propagation speed 
approaches vp in the direction normal to the field line 
(0 = 90°) and drops to the Alfven speed when = 0°. 
But, as a function of height, these two speeds become 
equal and the fast wave speed approaches the Alfven 
speed. The dark-grey shaded regions (Regions IV, V 
& VI) correspond to the propagation of slow waves on 
the axis of the flux tube. The slow wave has a maxi- 
mum propagation speed when = 0° and in this region 
it is equal to the sound speed. The propagation speed is 
slower than sound speed in any other direction (0 > 0°) 
and vanishes to zero in the direction normal to the field. 
Hence, in the case of the propagation along the axis of 
the flux tube, we notice a distinct separation between 
propagation of slow and fast waves and also that there 
is a region where there is no wave propagation possible. 



This scenario remains the same even as one moves away 
from the axis of the flux tube. 

Ambient medium: Let us now look at the propagation in 
the ambient medium, represented here by regions III and 
VI (shown in Fig. 6 by vertical-line pattern). Region III 
is the propagation zone of the fast wave in the ambient 
medium. The speed in the direction normal to the field 
line is vp which is close to sound speed in the lower part 
of the medium but increases with height as the Alfven 
speed increases. The propagation at = 0° is set by 
the sound speed below the height marked as B at which 
point the sound speed and the Alfven speed are equal 
and from there upwards the propagation in the direction 
of the field is set by the Alfven speed. The slow wave 
propagation in the ambient medium is marked as dark- 
shaded region with horizontal strips (Region VI). Up to 
the height £>, the Alfven speed (solid grey line at the 
bottom left part of the figure) is lower than the sound 
speed and hence the speed of propagation along the field 
line (the maximum) is the Alfven speed but above the 
height B it shifts to sound speed. In essence, we do not 
expect any wave propagation in the ambient medium cor- 
responding to regions IV and V. The fast and slow waves 
are distinctly separate below the height of point B. 
Along field line: Let us analyze the propagation speeds 
on a field line that crosses the height z = 1 Mm at radial 
distance of r = 0.35 Mm from the axis and azimuth of 
1.757T radians. The propagation region of the fast wave 
on this field line is marked in Fig. 6 as Region II and 
the slow modes are present in Regions V & VI. Unlike 
in the ambient medium, the fast wave propagation nor- 
mal to the magnetic field (equal to vp) is larger and we 
should start seeing an expanding wavefront with height 
above 0.3 Mm. The propagation speed in the direction 
of the field remains cs (thick dot-dashed curve) till point 
A from whereon the propagation speed is va (thick solid 
curve). The slow waves on this field line exist in regions 
V and VI (marked by grey-shade and slanted lines). Up 
to a height of point A, the maximum speed of propa- 
gation is Alfven speed but it changes to sound speed 
as the perturbations crosses this height. Similar to the 
above mentioned cases of propagation along the axis and 
ambient medium, there exists a region where we do not 
expect any wave propagation (Region IV). Due to the 
increasing Alfven speed with height, the speed of fast 
mode is approximately equal to the Alfven speed on the 
flux tube axis and hence distinguishing between the in- 
termediate Alfven mode and the fast mode is in practice 
rather difficult, unless polarization information is avail- 
able. In the ambient medium the difference between vp 
and va is slightly larger, making it feasible to separate 
out the three modes. The excitation mechanisms, that 
we consider in this paper, generates magneto-acoustic 
waves within the magnetic flux tube as a result of the 
driving motion at the footpoint. Since we do not have a 
variation within the magnetic flux tube from high-/? to 
low-/? plasma, we cannot separate the energy contribu- 
tion from the fast and Alfven modes. 

3. WAVE EXCITATION 

The study of the excitation of waves in a vertically 
oriented magnetic flux tube embedded in an otherwise 
unmagnetized ambient medium has been discussed in a 
series of papers by Musielak et al. (Musielak et al. 1989, 
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Fig. 6. — Different propagation regimes of the magneto-acoustic 
waves present in the flux tube as function of height. The plot 
shows the characteristic speed of waves on the axis of the flux tube 
(shaded grey: Region I, IV, V & VI), on the representative field 
line (slanted line patterns: Region I, V & VI) and in the ambient 
medium (horizontal line patterns: Region III Sz VI). The character- 
istic speeds as a function of height are also shown: Alfven (solid), 
fast (dashed), tube (dotted) and the sound speed (dot-dashed). 
The thick curves correspond to the values on the representative 
field line. 

1995, 2000, 2002). They have outlined the theory of the 
interaction of a flux tube with turbulent convection and 
the consequent generation of waves in the magnetic fea- 
ture for solar and in general for stellar atmospheres. Ex- 
citation by granular buffetting, a process by which waves 
can be generated in magnetic flux tubes have been stud- 
ied by Hasan et al. for a flux tube in a photospheric net- 
work (Hasan & Kalkofen 1999; Hasan et al. 2000), and 
for explaining the heating of the chromospheric network 
(Hasan et al. 2005; Hasan & van Ballegooijen 2008). 

In a way similar to Hasan et al. (2005) and others 
(Hasan & van Ballegooijen 2008; Vigeesh et al. 2009; Fe- 
dun et al. 2011b), we investigate wave propagation by a 
range of excitation in the equilibrium model. The excita- 
tion is driven at the lower boundary by prescribed drivers 
that model different wave generation mechanisms at the 
photospheric layer. Our earlier work (Vigeesh et al. 2009, 
2011) considered wave excitation in a 2D equilibrium 
magnetic field configuration by a transverse motion of the 
lower boundary similar to Hasan et al. (2005). The possi- 
ble wave modes generated in this scenario were restricted 
to the slow and fast magneto- acoustic modes. Since we 
now turn to a 3D configuration in this paper, studying 
the properties of yet another kind of wave, namely, the 
intermediate Alfven wave is now also possible (Fedun 
et al. 2011b). Nevertheless, it would be worthwhile to 
look into the 3D problem of granular buffetting, which is 
not likely to generate any Alfven waves. 

3.1. Horizontal driver 

We implement a horizontal driver to mimic the granu- 
lar buffetting motion as observed in photospheric layers. 
The transverse driving velocity at the bottom boundary 
is specified as follows: 

y,(x,„,o, t) -[**•{¥) °^ p / 2 ' 

[ t > P/2. 



The amplitude of the motion (Vb) is 750 m s" 1 and wave 
period (P) is 24 s which is quantitatively in the range of 
values found in subsurface convective layers in which the 
flux tube is rooted (Berger & Title 1996; Utz et al. 2010; 
Keys et al. 2011). 

3.2. Torsional driver 

In order to study the effect of vortex-like motion in a 
flux tube, we also use a torsional driver. In this case, 
the azimuthal component of the velocity is specified as 
follows: 

V,( X , y, o,t) = { ~ v ° tanh (?) sin (?) * ^ p ^ 

[0 t>P/2. 

Here Sr describes the characteristic distance by which 
the azimuthal component changes from to Vq. For a 
practical purpose we use Sr = 1 Mm. The amplitude 
and period in this case are the same as for the horizontal 
driver. 

The driving motion prescribed above generates a wide 
range of wave-like disturbances in the medium which we 
analyze in the following section. 

4. NUMERICAL METHODS AND BOUNDARY 
CONDITIONS 

The 3D numerical simulations were carried out using 
the Sheffield Advanced Code (SAC) (Shelyag et al. 2008). 
SAC is a fully non-linear MHD code based on a more 
general-purpose code named VAC (Toth 1996). SAC 
uses a modified version of the set of MHD equations to 
deal with a strongly-stratified magnetized medium. The 
spatial derivatives are calculated using a fourth-order 
central difference scheme and the time derivatives are cal- 
culated using a Runge-Kutta scheme. Various boundary 
conditions are available for the different physical prob- 
lems. For our purpose we use open boundary conditions 
for the top and the side boundaries. This allows any gen- 
erated perturbation to propagate out of the simulation 
domain without reflection. 

5. SIMULATION RESULTS 

Here we have considered a strong magnetic flux tube 
with the plasma- f3=l surface outlining the boundary of 
the constructed magnetic flux tube. Since, given that, in- 
side the magnetic flux tube plasma /3 < 1, the magneto- 
acoustic waves generated by the driver are the slow (pre- 
dominantly acoustic) and fast (predominantly magnetic) 
waves (see Vigeesh et al. 2009). 

The gas pressure perturbations are directed along the 
field lines. The slow magneto-acoustic wave (SMAW) in 
a low-/3 plasma can be visualized by the parallel compo- 
nent of velocity (vu) to the background magnetic field, 
as they are predominantly acoustic in nature. On the 
other hand, as the fast magneto-acoustic wave (FMAW) 
can also propagate across the magnetic field lines, it can 
be identified in the perpendicular components of veloc- 
ity (v±) to the background magnetic lines of force. The 
propagation of the intermediate mode is characterized by 
the azimuthal component (v^) of the velocity. The de- 
scription of the three components of velocity is shown in 
Figure 7. 
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In the following, we describe the method used to 
compute the three components of velocity. Since the 
initial model is axisymmetric, we can assume that the 
cross-section of the iso-surface, or the surface of equal 
magnetic potential of the flux tube at a given height 
traces a circle with the center at axis. See Fig. 2 for 
more information on the geometrical shape of the flux 
tube and the position of the field line. We proceed by 
calculating the field lines that pass through these points 
at time t = 0. The displacement of the starting points 
are determined at each timestep and the new magnetic 
field lines are traced as they evolve. The driving motion 
at the base distorts the flux tube and consequently the 
flux tube departs from being axisymmetric. However, 
the field lines that define the magnetic flux iso-surface 
remain the same. For brevity, we consider a single field 
line at a radius of r=0.35 Mm and azimuthal angle of 
<p = 1.757T radians (dashed curve in Fig. 2 & Fig. 7). 
The parallel component of velocity (vn) along the field 
line (s) can be computed by taking the scalar product of 
v on B. We construct the vector normal to the surface 
of equal magnetic potential (n) by using two adjacent 
field lines on the points that form the circle. The set of 
points which define the two adjacent field lines form a 
plane and we use these points to calculate the surface 
normal, h. The projection of the velocity vector (v) 
on this vector gives the normal component. The cross 
product of the magnetic field vector (B) and the normal 
vector then gives the vector tangent to the magnetic 
equipotential surface ((/>). The projection of the velocity 
(v) on this vector gives the azimuthal component, v^. 

The parallel (s), normal (n) and the azimuthal (</>) 
vectors for a point P on the selected field line (thick 
dashed line) are shown in Fig. 7. The magnetic field 
lines that form the iso-surface are shown as set of thin 
lines. Determining the velocity components will help in 
separating out the contribution of different modes that 
are generated as a result of the driving motion at the 
lower boundary. 

In the following we discuss the result of our numerical 
simulation for the two drivers discussed in the previous 
section. 

5.1. Horizontal excitation 

We consider the case in which the magnetic flux tube 
is excited by a uni-directional horizontal motion at the 
lower boundary. The physical scenario is similar to the 
case wherein a deep rooted flux tube in the interganular 
lane is buffeted by surrounding granules. From the dy- 
namics of the magnetic bright points (MBPs) , as studied 
by Cranmer & van Ballegooijen (2005), this motion cor- 
responds to a single instance wherein the magnetic flux- 
tube receives a random kick from a direction not neces- 
sarily the same along which it was already moving. Such 
a transversal driver generates strong SMAW and FMAW 
in the medium that propagate along the magnetic flux 
tube. Figure 8 show the time-distance plots of the ve- 
locity component as described in the beginning of this 
section. These panels correspond to vn , v± and along 
a vertical path parallel to and about 140 km away from 
the magnetic flux tube axis. 

The longitudinal perturbations produced by the pas- 
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Fig. 7. — The description of the parallel (v\\), perpendicular (v±) 
and the azimuthal (y^) components of velocity vectors with respect 
to the magnetic field line and the equipotential surface. 

sage of the slow wave along the field line can be visualized 
by plotting the parallel component of velocity (v\\) on 
the flux tube magnetic field. We find dominant velocity 
perturbations associated with the SMAW. This driver 
also generates relevant FMAW in the medium. Since 
the FMAW can also propagate across the magnetic field 
lines, we determine the perpendicular component of ve- 
locity (v±) to the magnetic field of the flux tube. The 
azimuthal component of velocity (v<j>) has a comparable 
magnitude, as seen from the contour values. This is a 
result of the location of the particular field line under 
consideration with respect to the driving motion which 
results in a significant contribution to the azimuthal com- 
ponent. 

Figure 9 shows a snapshot (at t=100 s) of the tem- 
perature (left panel) and the magnetic field (right panel) 
fluctuations on x — y planes at three different heights 
(z=150 km, 500 km and 800 km). Few field lines mark 
the magnetic flux tube and the variation of the field 
strength along the field lines are color coded. The snap- 
shot depicts a clear signature of the passage of the SMAW 
at height z =500 km. Glancing through the sequence of 
the snapshots, in temporal order: 

• at t~35 s, strong velocity and the magnetic field 
perturbations suggest the arrival of the fast MHD 
mode at z =500 km; 

• at t~46 s, the fast MHD mode crosses at z =800 km 
level resulting in amplified velocities and magnetic 
field perturbations. The circular contours corre- 
sponding to a constant magnetic field strength val- 
ues, shown on top of the box, shifted with respect 
to the time at which the FMAW reaches the top of 
the box; 

• at t~67 s, the slow MHD mode crosses at 
z =500 km causing strong temperature fluctuations 
(ST); 

• similarly at, t =100 s (Fig. 9), there is a strong 
evidence of ST which is absent in the case of a 
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Fig. 8. — The time distance plot of the three components of velocity (v\\,v_\_ and v^) along the vertical path slightly away from the axis 
of the magnetic flux tube as a result of horizontal excitation. 



torsional excitation. 

Since the flux tube has a substantial B x component at 
the footpoint as seen in Fig. 1, the velocity of the uni- 
directional driver has a larger and more dominant com- 
ponent parallel to the field in a small region on either side 
of the magnetic flux tube. We see strong v\\ components 
on either side of the magnetic flux tube concentrated on 
the axis associated with the SMAW. 

5.2. Torsional excitation 

Observations of the solar photosphere reveal that 
bright points corresponding to flux tubes located in the 
intergranular lane may follow a spiral path centered at 
regions where boundaries of two or more granules meet 
(Bonet et al. 2008). One assumes that an effective torque 
acts on these flux tubes, and they start to rotate as they 
follow such a trajectory. The azimuthal component be- 
comes amplified, as these flux tubes approach the strong 
down-flow regions at the confluence of two or three in- 
tergranular lanes. Assuming that the simulated mag- 
netic flux tube is located where the azimuthal driving is 
stronger, we study the dominant wave modes that are 
likely to be excited. 

Figure 10 shows the time-distance plot of v\\ : v± and 
components of velocity where the MHD perturbations 
are driven by a torsional driver. The vn amplitude in the 
case of torsional excitation is rather small (~30 ms" 1 ) 
compared to the horizontal excitation (~300 ms" 1 ), sug- 
gesting that the torsional motion generates SMAW that 
are weaker than those generated by a horizontal driver. 
This feature is also seen in the temperature fluctua- 
tions, which are mainly associated with the slow MHD 
mode. Also, the torsional driver generates weaker tem- 
perature fluctuations than the horizontal driver. The 
time-distance plot of vn shows a weak fast MHD mode 
branch that splits from the dominant slow mode branch 
which is absent in the case of a horizontal excitation. 

We noted above, in the case of horizontal excitation 
(Sec. 5.1), that the actual values are calculated on a verti- 
cal line parallel to the magnetic flux tube axis. Since the 
magnetic flux tube expands with height, the represen- 
tative vertical line crosses a range of different magnetic 
equipotential surfaces, significantly in the lower part of 
the computational box. The propagation of the interme- 



diate mode is characterized by the rotation of magnetic 
equipotential surfaces with different amplitudes as one 
moves away from the flux tube axis. We have found that 
the magnetic flux tube suffers a torsional kick and this 
motion is carried away along the field lines. The contours 
of constant magnetic flux value rotate independently of 
each other as can be seen at the top of the box. Even 
though the B x component at the foot point of the mag- 
netic flux tube is significant, the velocity of the torsional 
driver is in the plane normal to the magnetic flux tube 
and hence nowhere parallel to the field. Thus, SMAW 
are not excited directly and dominantly in this case. 

5.3. Energy transport 

Let us now consider the energy transport by MHD 
waves for the two driving cases detailed in the previous 
section. Following Vigeesh et al. (2009), instead of using 
the full non-linear equation for energy flux, we consider 
the linearized wave-energy flux equation given by Bray 
& Loughhead (1974). The acoustic (Fa) and Poynting 
fluxes (Fp) are computed for the two cases as follows: 

F A = ApV, (3) 



F P = -?-(B -AB)V- 

47T 



47T 



(V-AB)B , 



(4) 



where A represents the perturbations with respect to 
the equilibrium and Bo is the equilibrium magnetic field. 
Comparing the acoustic flux generated by the horizontal 
and torsional drivers, we find that in the case of a hor- 
izontal driver, there is a significant amount of acoustic 
flux crossing the height z =500 km. We already discussed 
in Sec. 5.1, the uni-directional driver generates strong 
SMAW. Hence, we see a strong acoustic flux (longitudi- 
nally dominated) associated with it. This is absent in 
the case of torsional driver, as there is negligible amount 
of slow wave generation (see Sec. 5.2). 

Since the slow, fast and torsional Alfven modes are 
present in the medium, it is inspiring to calculate the 
contribution of these three modes in terms of the wave 
energy that they transport. The fast and slow MHD 
mode contribute to the acoustic (Fa) and Poynting fluxes 
(Fp) over different scales depending on the plasma- /?. 
Assuming that the driver implemented here generates 
linear waves, we have calculated the wave energy fluxes 
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Fig. 9. — Temperature (left) and magnetic field strength (right) fluctuations in a flux tube at three different heights [z =150 km, 500 km 
and 800 km) at t=100 s as result of a transverse uni-directional excitation. Field lines are shown to mark the flux tube. The contour of 
equal magnetic field strength at z =800 km are shown on the top. 
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Fig. 10. — The time distance plot of the three components of velocity (t>||,i>_i_ and v$) slightly away from the axis of the magnetic flux 
tube as a result of torsional excitation. 



of the different modes for two cases. For any given mode, 
the energy of the wave motion transported through the 
medium depends on the direction with respect to the 
magnetic field and maximum speed of the mode in the 
particular direction. Decomposing these energy compo- 
nents in the three directions, the wave energy fluxes can 
be computed as follows, 



1 = P v \\ c s, 



pv\\lc\ 



(5) 



(6) 



(7) 



where, Fn , F± and F^ are the parallel, perpendicular and 
azimuthal energy fluxes corresponding to the low plasma- 
/3 propagation of slow, fast and intermediate wave, re- 
spectively. It should be noted that the perpendicular 
component of the time-averaged energy computed by 



Khomenko & Cally (2012) is slightly different from the 
one that is used here. The maximum speed of the FMAW 
in the perpendicular direction to the magnetic field is 
equal to vp. Although vp is equal to the Alfven speed 
at higher layers on the axis of the flux tube, this is not 
the case for the particular field line that we have chosen 
as shown in Fig. 6. Hence, it is appropriate to use vp to 
calculate the perpendicular wave energy flux. 

Figures 11 and 12 show the time-distance plot of the 
three energy fluxes (Fn, F± and F^) computed on the 
representative magnetic field line as a result of horizon- 
tal and torsional excitation in a flux tube, respectively. 
At larger heights, the parallel flux (Fy) corresponding to 
predominantly SMAW for the case of uni-directionally 
excited flux tube is O(10 7 ) erg cm -2 s _1 , which is two 
orders of magnitude larger than for the case of a torsion- 
ally excited flux tube. There is a significant amount of 
perpendicular flux, (9(10 6 ) erg cm -2 s _1 , being trans- 
ported close to the Alfven speed in the case of horizontal 
excitation which is also absent in torsionally excited flux 
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tube. As noted in Sec. 5.1, the contribution to the az- 
imuthal component as a result of the location of the flux 
tube with respect to the driving motion results in a sig- 
nificant azimuthal energy flux for the unidirectional case. 
There is indeed a more stronger azimuthal energy flux for 
the torsionally excited flux tube. 

6. SUMMARY AND CONCLUSION 

Three-dimensional modelling of the MHD wave prop- 
agation in a magnetic flux tube embedded in a highly 
stratified solar atmosphere was carried out. We have in- 
vestigated two types of excitation mechanisms that are 
characteristic processes by which waves can be generated 
in magnetic flux concentrations of the real solar atmo- 
sphere. We considered a transversal driver that mimics 
granular buffeting and a torsional driver which can be ap- 
proximated to observed vortex flows within flux tubes. 
The driving motion generates slow and fast magneto- 
acoustic waves in the medium. The granular buffeting 
seems to generate strong SMAWs and FMAWs, when 
compared to vortex-like motion. The longitudinal acous- 
tic flux transported by these waves also varies depending 
on the excitation mechanism. We found that, as a result 
of the granular buffetting simulated here using a unidi- 
rectional horizontal motion of the lower boundary, more 
acoustic energy is generated and this is likely to become 
dissipated at chromospheric heights. Hence, we conclude 
that the granular buffeting produces stronger emission in 
the chromospheric network when compared to the vortex 
flows driven by convective downflow. The observations 
of vortex flows in the photosphere and the corresponding 
chromospheric feature reveal that the bright points asso- 
ciated with the flux concentration follow a spiral path to- 
wards the center of the vortex. Although high-resolution 
observations do not completely resolve these features yet, 
one may realistically suspect that these bright points are 
constantly buffetted by the granules on their way to the 



intergranular sink. We suggest that these buffeting re- 
sults in a horizontal excitation of the magnetic flux tube, 
generating slow longitudinally dominant MHD waves and 
FMAW which eventually dissipate due to shock forma- 
tion and show up as emission in the chromosphere. How- 
ever, as these flux tubes approach the eye of the vor- 
tex, the horizontal excitation ceases and the flux tube 
is mainly driven by a torsional motion, which generates 
negligible SMAWs and FMAWs. This results in lesser 
shocks in the upper atmosphere and, hence, a weaker 
emission. 

Apart from the other sources of wave driven heat- 
ing mechanisms like, strong downflows in the outskirts 
of magnetic elements resulting in excitation of slow 
modes within the flux tube developing into shocks (Kato 
et al. 2011), transversely shifting magnetic field lines 
within the magnetic concentration producing Alfven 
waves causing Alfven wave turbulence in the upper atmo- 
sphere (van Ballegooijen et al. 2011), the simple process 
outlined in this paper is suggested to be considered as a 
viable mechanism that can contribute to the heating of 
magnetic elements. To be more precise, we conjecture 
that more than one mechanism of wave production is ef- 
fective during the excursion of a magnetic element in the 
erratic intergranular lane, and result in temporal varia- 
tions of the overall emission that occur at chromospheric 
heights over these magnetic elements. 
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483133). We would like to thank the anonymous referee for 
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APPENDIX 

INITIAL MAGNETO-HYDROSTATIC MODEL 

In this section we develop a 3D model for the magnetic field B(r) and plasma parameters of a photospheric 
magnetic flux tube. A Cartesian reference frame (x, y, z) is used with z the height above the base of the photosphere 
(i.e., the level where the optical depth T5000 = !)• The computational domain is given by x = [— y = [— 
and z = [0, L z ], where L and L z are the half-width and height of the domain, respectively. The axis of the flux tube 
is vertical and lies along the z-axis of the reference frame. The flux tube is symmetric with respect to the planes 
x = and y = 0, so for the purpose of computing the initial equilibrium field only one quarter of the domain needs 
to be considered: x = [0,L] and y = [0,L]. At low heights the flux tube has a nearly circular cross-section with 
radius tq < L, and is surrounded by a nearly field-free medium. However, the gas pressure decreases with height z, 
causing the cross-section of the tube to increase with height, and at larger heights the radius of the flux tube becomes 
comparable to the half-width L of the domain. To simulate the effect of the neighboring flux tubes, we assume that 
the normal component of B vanishes at the side walls of the domain, so that the field lines are forced up along 
the side walls. In this paper we use L = 0.5 Mm, corresponding to a total width of the flux tube of 1 Mm at large height. 

The magnetic field B(r) is written in terms of Euler potentials: 

B = Vu x Vu, (Al) 

where u(r) and v(r) are constant along field lines (B • Vix = B • = 0). At the base of the model the vertical 
component of the field has a Gaussian profile, 

B z (x,y,0) =B exv[-(x 2 + y 2 )/r 2 ], (A2) 

where ro is the 1/e- width of the profile (we use ro = 120 km) and Bq is the peak field strength. The Euler potentials 
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Fig. 11. — The time distance plot of the three energy fluxes (Fy ,Fj_ and F^), corresponding to the slow, fast and intermediate wave for 
a distance along the vertical line slightly away from the axis of the magnetic flux tube as a result of horizontal excitation. 
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Fig. 12. — The time distance plot of the three energy fluxes (Fy ,F± and F^), corresponding to the slow, fast and intermediate wave for 
a distance along the vertical line parallel to the axis of the magnetic flux tube as a result of torsional excitation. 



at the base are chosen such that u(x, y, 0) depends only on x and v(x, y, 0) depends only on y: 

u(x, y, 0) = G{x) and v(x, y, 0) = G(y), 

where 

G(x) = x/Bq I exp ( — — 7T ) ds, 
Jo V r$J 

which is consistent with equations (Al) and (A2). The side boundary conditions are 

u(0, y, z) = 0, u(L, y, z) = G max , v(x, 0, z) = 0, v(x, L, z) = G max , 



(A3) 
(A4) 

(A5) 



where G max = G(L). At the upper boundary no conditions on u(x,y,L z ) and v(x,y,L z ) are imposed, allowing the 
field lines to freely move on the upper boundary. Hence, there should be no build-up of field-aligned electric currents 
inside the domain. The interior of the domain is assumed to be in magnetostatic equilibrium: 



+ pg - 



1 

47T 



(V x B) x B = 0, 



(A6) 



where g = — gz is the gravitational acceleration, p(r) is the plasma pressure, and p(r) is the mass density. The third 
term on the left-hand side describes the Lorentz force due to distributed electric currents at the interface between the 
flux tube and its surroundings. The gas pressure and density are given by 



p(u, v, z) =Pint(z) [1 + A) VO, v)] , 
p(u,v,z)=p int (z) [1 + /3 "V(^)] , 



(AT) 
(A8) 



where Pmt(^) and p[ n t(z) are the internal pressure and density along the axis of the flux tube, /3q is the ratio of gas- 
and magnetic pressures at the origin (z = 0), and f(u, v) varies smoothly from zero on the flux tube axis to f(u, v) = 1 
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in the external medium and at the side walls. The internal density is given by p[ n t(z) = —g~ 1 dpi nt /dz, so there is 
hydrostatic equilibrium along field lines. Note that p/p depends only on height z, so the temperature is constant in 
horizontal planes. The internal pressure as function of height is approximated as a sum of two exponentials: 

Pint (2) =Piexp(-z/H 1 ) +p 2 exp(-z/# 2 ), (A9) 

with a photospheric pressure scale height Hi = 110 km (pi = 4.2 x 10 4 dyne cm -2 ), and a chromospheric scale height 
H2 = 220 km (p2 = 10 3 dyne cm -2 ). The maximum of the magnetic field strength at the foot point is estimated using 
the thin tube approximation: Bo = y/%irpint{fy/ Po- To obtain kilogauss fields in the photosphere we use /3q = 0.5, so 
that the external gas pressure is three times the internal gas pressure. The function f(u,v) is given by 

f(u,v) = F(A) = (4A-A A )/3, (A10) 

where 

A(u,v) = l-cos( -^—) cos (-^-) . (All) 

Note that A(u, v) varies from A = on the axis (u = v = 0) to A = 1 in the field-free medium and on the side 
boundaries (u = G max or v = G max ), and that dF/dA = at A = 1. Hence, there is a smooth transition of the gas 
pressure from the interior to the exterior of the flux tube, and the electric currents at the interface are finite. 



Inserting expressions (Al), (A7) and (A8) into equation (A6) yields two coupled equations for the Euler potentials 
u and v. 

^(VxB).V,- Pint (^ - 1 g=0, 

^ ( V x B) • Vu + p int (z)^ 1 1£ = 0. 
The solution of these equations can be obtained by solving a variational problem for the Lagrangian, 



W 



rL r>L rL z 

Jo Jo Jo 



\Vu x Vv\ : 

8^ 



"Pint (*) A) 1 f( u > v ) 



dx dy dz, 



(A12) 
(A13) 

(A14) 



where the variations in (u,v) are subject to the above-mentioned boundary conditions. To solve this variational 
problem, we set up a grid of 50 x 50 x 100 cubic cells covering the interior of the computational domain, and ghost 
cells are added at the side boundaries. The Euler potentials and are defined at the cell corners (i, j and k 
are indices on the grid), and the derivatives in equation (A14) are evaluated using second-order finite differences. The 
Lagrangian W(uijk, vijk) is then a function of the potentials at all interior grid points. We minimize this quantity 
using the conjugent-gradient method (Press et al. 1992). Once the solution is obtained, the magnetic field, pressure 
and density are evaluated at cell centers, and the computed values are mirrored to obtain a 3D magnetostatic model 
of the entire flux tube. This model is used as the initial condition for the 3D MHD calculations. 
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